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The tomographic map of quantum state of a system with several degrees of freedom which depends 
on one random variable, analogous to center of mass considered in rotated and scaled reference frame 
in the phase space, is constructed. Time evolution equation of the tomogram for this map is given in 
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in details. Symmetry properties in respect to identical particles permutations are discussed in the 
framework of proposed tomography scheme. 
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I. INTRODUCTION 

Recently a new formulation of quantum mechan- 
ics was suggested. This new formulation uses non- 
negatively defined probability distribution function to de- 
scribe quantum states, called marginal distribution HQ, 
or tomogram. This function can be considered as an ana- 
log of known quasidistribution functions like nonnegative 
Husimi Q-function Q or Sudarshan-Glauber P-function 
H 0| • Tomographic approach was initially developed for 
one-mode systems, in this case the quantum state is de- 
scribed by the density matrix p(q',q") HQ or by sym- 
plectic tomogram w(X,fi,u) |{j|. Here density matrix is 
the function of 2 variables and tomogram is the function 
of 3 variables. Seeming overcompleteness of the tomo- 
graphic description is balanced by the fact that quantum 
tomogram is a homogeneous function [Tol 1 1 1| . 

The state of the system with N degrees of freedom is 
described by density matrix p(<f , if'), the function of 2N 
variables. Then, what is the generalization of quantum 
tomogram, will it depend on 3iV variables or on 2N + 1 
variables? Such generalization was developed for the to- 
mogram, depending on 3iV variables (usual symplectic 
tomography) |l2L Il3| - In this paper we propose the tomo- 
graphic map with only one random variable (i.e. 2N + 1 
variables totally) and discuss its properties in details. 

Quantum tomography becomes popular nowadays for 
a number of reasons. It was developed not only for con- 
tinuous variables like position but also for discrete spin 
variables [3 13 13 111 13 13 113 First, tomographic 
representation operates with the values, which are di- 
rectly measured in experiments, for example non-classical 
and coherent states of fehtJHM 13 111 13 13 111 113 
or matter optics [3 |3Cj, l3ll l32| experiments. Second, 
the tomogram is non-negative, and this must attract the 
attention of those who deal with computer simulations of 



quantum systems |33j . Many problems in this field, e.g. 
for Path Integral Monte Carlo simulations, result from 
the use of the alternating-sign, or even complex, values 
to describe the quantum state (for example, sign prob- 
lem in Fermi-systems simulations). The latter leads to 
difficulties with convergence which can be overcome by 
quantum tomography approach. These facts along with 
the tomography applications in quantum computations 
and entanglement (see, e.g., |34| ) . as well as in the the- 
ory of information and signal analysis |2ri| , show that the 
development of convenient and simple tomographic map 
for the case of many particles (or, which is the same, 
many modes) is a task of great significance. 

The paper is organized as follows. In Sec. [H] we present 
the definition of tomographic scheme with one random 
variable, elucidate some of its useful properties and dis- 
cuss the physical meaning of the map. In Sec. II I II wc 
derive the equations describing quantum evolution, sta- 
tionary states, quantum transitions and rules for aver- 
age values calculation for the proposed tomography map. 
Some examples of state description using the developed 
approach are given in Sec. II VI and symmetry of the map 
with respect to particles permutations is discussed in 
Sec. The work is summarized in Sec. I VII 
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II. ONE RANDOM VARIABLE TOMOGRAPHY 

A. Definition of the tomographic map 

We begin with the one-dimensional (ID) case of a par- 
ticle with continuous degree of freedom (in this paper 
we do not consider spin variables, but generalization of 
the formalism is straightforward). Quantum mechanics 
states that we know 'everything' about the system if we 
know density matrix. In practice, to obtain any informa- 
tion about the system we have to measure some quanti- 
ties, for example coordinate q or momentum p. It is also 
possible sometimes to measure an intermediate quantity, 
fiq + vp, where /i, v are real parameters. Formally, this 
quantity (denote it X) is coordinate, measured in scaled 
and rotated reference frame in the phase space. It turns 
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out, that the distribution function of X (w(X, p,, i/)), 
measured for all sets of jj,, v gives complete quantum me- 
chanical description of the system, in the sense that there 
is a unique correspondence between w(X,p, v) and den- 
sity matrix (see, e.g, 0, |U El EE E3 ) Note that distri- 
bution function w(X, p,, v) is equal to (5(X — pq — vp)}, 
where (...) is quantum mechanical average. Then there is, 
in principle, a possibility of complete experimental den- 
sity matrix determination through the set of coordinate 
measurements. 

When we deal with more than one particle and di- 
mension we can consider individual Xj = pjqj + vjpj 
for every j-th degree of freedom. This results in the 
symplectic tomography representation 0, . Here we 
are to show that it is enough to work with only one 
X = J^jXj- To do this, let us consider the general- 
ization of w(X, p, v) = (8(X — pq — vp)) i where q,p and 
p, v becomes the vectors, their products become scalar 
products of vectors, while X remains a real number: 



D(y) = {2ir)- Nd exp 



w(X, ft, v) = (S(X — jtq — vp) 



(1) 



Related problems were discussed in (271, |28| . 

Throughout the paper designations are the following. 
We consider the system of N particles in d dimensions, 
the number of degrees of freedom is Nd. Vectors are writ- 
ten as a, we use everywhere the vectors with Nd com- 
ponents, if the otherwise is not mentioned. Designation 
e is used for the vector with all components equal to 1 
(cj = 1). Scalar product of vectors is designated a = be 



(meaning a = ^ . bjCj 



hoc denotes the compo- 



nentwise product of vectors (aj = bjCj). The tomogram 
for usual symplectic scheme is designated as w s (X, ft, v) 
{X,jl and v with Nd components each); the tomogram 
with one random variable is written as w(X, ft, if) ■ We 
also use Planck constant h = 1 everywhere. 

We can begin the construction of one-random-variable 
tomography representation from the known star product 
expressions. In the framework of star product formalism 
HE EM 03 EH every operator is replaced by the function 
('symbol'), depending on specific set of parameters (y), 
products of operators turn into 'star products' of corre- 
sponding symbols (star product problem was discussed 
also in Ref . |42l| ) . As a result, one deals with functions 
only, avoiding operators. For example, using a pair of op- 
erators D(y), U(y), we construct the connection between 
the symbols f A {y) and operators A: 



f A {y) = Tr(AU(y)), 
A I f A (y)D(y)dy, 



(2) 
(3) 

Tr(D(y)U(y))dy = 1 (4) 
For y = {X, ft, v} one can choose 

U(y) =S(X- M- vp) , (5) 



i ^X — flq '— vpj 



(6) 



which defines the symbols (denote them wa{X, ft, v)) and 
star product: 

(w A *w B ){y) = WA(y")w B (y')K(y",y',y)dy"dy' (7) 



The kernel of star product K(y" , y', y) is expressed as 
follows: 

K(y",y',y)=Tr[D(y")D(y')U(y)] = 

e - i(k x-x'-x") § ^ + j _ h ^ b{ p, + if _ k0) x 

i{Ji" v — k{jt" -\-p? + v' v' 1 +k 2 £Lv) /2) 



(27T) 

For any operator A we have (A) = Tr(pA), therefore 
Wp-symbol of density operator p is the same as w(X, ft, v) 
defined by Eq.Q. Density matrix in any representation 
is just a matrix element of density operator. Then, fi- 
nally, we come to the unique correspondence (invertable 
map) between the tomogram w and density matrix men- 
tioned above: 



(9) 



w(X, ft,v) = (5(X - flq- vp) ) , 



r , w {X,ft,v)e^-^^^- (10) 



(27T) 



Density matrix always can be reconstructed from the to- 
mogram w using these equations, so one random variable 
tomogram describes quantum state completely. Note that 
now the state-describing function is nonnegative and de- 
pends on 2Nd + 1 variables in contrast to symplectic 
tomogram, depending on 37V d variables. 



B. Properties of w(X,j2,v) 

It is convenient to consider density matrix in coor- 
dinate representation and Wigner formulation of quan- 
tum mechanics ,35] to derive the properties and evolu- 
tion equation for the tomogram w. In the framework of 
Wigner formalism the state of the system is described by 
real Wigner function W(q,p) defined in phase space and 
connected with the density matrix as follows: 

nm = J P(r+ f f^^p- (") 

p(g',q") = J W (^±^,p)e^-^dp (12) 
Using Eas. lPTU|) and Eas. (|llll2[> we obtain: 

w{X,ft,v)= I W(q,p)e- ik ^-^-^^^,(13) 



W(q,p) 



e -Hm-0P-x) w (x, ft,v) 



(27T) 

dXdfldv 
(2ir) 2Nd 



(14) 
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Usual symplectic tomography map is developed in 
references [T3 and [3(1 . The symplectic tomogram 
w 8 (X,jl,u) and Wigner function are connected as fol- 
lows: 



w s (X,p,u) 



W(q,p) 



W(q,p)e 



-ik(x-jloq-vop) dkdqdp 



(2tt)' 



-^^-^w^P^^SS (16) 



(27T) 



Since the Wigner function is connected by invertablc 
maps with both tomograms w and w s it is obvious that 
they contain the same information about the quantum 
state. In fact one has 

. Nd 

w{X,p,v)= w s (Y,p,v)5(X-J2Yj)dY, (17) 
J i=i 

w s (X,p,u) = J w{Y,kofi,kov)e t{Y - ijl) dkdY (18) 
The Wigner function is normalized: 
W(q,p)dqdp - 



,_ u _ u. .~ 7 dudqdp 



)5{u)dudq — I p(q,q)dq — 1(19) 



where we choose the normalization for density matrix 
Tr(p) = 1. Then the tomogram w is normalized in X 
variable: 



w(X,p,P)dX= / W(q,p)S(k)e lHM+pf > dkdqdp = (120) 



Although the tomogram depends on 2Nd+ 1 variables, 
instead of 2Nd for density matrix, the completeness of 
physical description is the same for both formulations, 
due to the fact that the tomogram is a homogeneous 
function. Consider the definition 113(1 and multiply all 
variables in w by a real number A: 



w{\X, A/t, Xv) = 

-ik(x-gff-aff) dkdqdp 



W{q,p)e- lXk{x - M -tp) dkd $ d P 



W(q,p)t 



(2tt|A|) 



(27T) 

w(X, /t, v) 



(21) 



where we just made the change of variables Xk — > k. 
Property l|21l) make obvious the relations 

w(X, p, V) = \X\~ x w(l, p/X, v/X) (22) 

For the pure state with wave function ^{q) symplectic 
tomogram w s was expressed in terms of modulus squared 
of fractional Fourier transform of the wave function (4^ • 
The tomogram w for pure state is given by 



w{X,p,v) 
^(q)exp 



dY 



S{X-^Yj] 



_X qo qY 



dq 



(23) 



C. Physical meaning 

We have defined the nonnegative function w(X,p,v) 
© completely describing quantum state. For any set 
of {p, v} it is normalized as a function of X, therefore 
w(X, p, V) is the set of distribution functions of quantity 
X. Then to know the quantum state completely one has 
to consider all sets of {p, u} (in practice, moving with 
some step) and measure X — pq + vp many times for 
each set: this yields the distribution function w(X, p,v) 
for given set of {p, z?}. 

Looking at Eq. l|22|l we see that we even do not have 
to know w(X, p,u), the value of this function in some 
point in X for all {p, v} is enough. This does not change 
the scheme of measurements, we still need to measure full 
distribution function of X for given {p, v\ (it is necessary 
to compare the values of distribution function in different 
points to be sure that statistical precision is good), but 
one has to store the smaller arrays of information. 

Property (|21|l can be used in another way. If A is equal 
to pp+Vv, we can parameterize {p, v\ by A and 2Nd— 1 
angles (to use the spherical coordinates in the space of 
{P,u}). Applying Ea. l(2T)l we come to reduced tomo- 
gram with 2Nd variables and {p, v\ located on the sphere 
with radius equal unity in 2A^<i-dimensional space. This 
new tomogram also completely describes the state and in 
some cases it can be convenient to use this one in mea- 
surements, because it is easier to sample 2Nd — 1 angle 
than 2A^d-dimensional space from —00 to 00 (see, e.g., 
[26l |28|). On the other hand such formulation causes 
trouble with the derivation of evolution equations and 
arbitrary average values calculation. 

The only remaining unclear point is the meaning of 
X = pq + up. It is the sum of positions measured in 
scaled and rotated reference frame in the phase space. 
But what does it mean physically? It is impossible to 
measure q and p simultaneously, but sometimes one can 
transform q and p into the form pq+vp, for example mix- 
ing the signal beam with local oscillator field (in quantum 
optics, see |25| and references therein). Another scheme 
was proposed in |26| , where q and p are mixed due to 
wave (electromagnetic or matter) propagation through a 
lens (or an analog of a lens in atomic optics). Taking 
into account the present development of science concern- 
ing controlling the Bose-condensates of atoms, this also 
can be a possible realm of tomography measurements. 
Bose-condensate is a coherent macroscopic state of many 
atoms and it is described by macroscopic wave function. 
For example, one can mix two such waves (condensates 
of the same atoms), using the first as a signal wave and 
the second as local oscillator. Varying the phase differ- 
ence of the condensates we sample different p, v. Prob- 
ably, the same can be done in superconductors (where 
the electrons of superconductivity also form the coherent 
macroscopic matter wave), using Josephson junctions. 

If we somehow accomplished the scaling and rotation 
of reference frame in the phase space we can measure 
the set of positions in this reference frame Xj = fijqj + 
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i>jPj , but it is enough to measure their sum, X = jlq + 
vp. It is analogous to the the position of center of mass 
measurement (the sum of coordinates of corresponding 
vector, to be more precise). Indeed, the center of mass 
position is 

X cm = 22 nijXj/M = ^ m.U'i'l, + VjPj)/ M , (24) 



where M — sumjrrij and rrij is the mass corresponding 
to j-th degree of freedom, and X cm can be associated 
with X = jlq + vp for some other set of {ft, v}. We 
sample all sets of {/I, v} therefore it is enough to measure 
the center-of-mass position in each scaled and rotated 
reference frame. 

Finally, we would like to make the following remark. 
The storage of arrays representing full density matrix or 
tomogram becomes impossible when the number of de- 
grees of freedom growth. If we use some grid, the num- 
ber of arrays elements is proportional to n Nd , where n is 
the number of grid steps. Increasing Nd we soon come 
to the situation when all data carriers in the world can 
not store corresponding arrays. And this is not neces- 
sary as the state of the system is uniquely determined by 
the one-particle density (through the density functional, 
see |44j and references therein). Then for many-particles 
systems description we can use reduced density matri- 
ces (one-body, two-body, etc.), and tomography map is 
constructed for them in the same way as for full density 
matrix. Then the situation with reference frame scaling 
and rotation is simplified because /j, and v are the same 
for all particles (if one-body density matrix is considered) 
and distribution functions are averaged over all particles. 



III. STATE TRANSFORMATIONS 



where p/rh means the vector with components Pi/rrii (the 
equation holds for the case of different masses for different 
particles and directions), the operators in the potential 
V designates the analytical expansion of the potential 
and use of the products of corresponding operators. This 
equation can be easily obtained applying the transform 
(JTTJ to the Eq.(|2J). 

To derive the evolution equation for tomogram one ap- 
plies the transform H13|l to evolution equation for Wigner 
function l|26|) . Expanding the potential in Ea. l26|) . it is 
seen that we have to consider the transforms of the fol- 
lowing quantities: qW, dW/dq, pW and dW/dp. The 
transform l|13fl of qW is 



qW(q,p)exp[—ik(X — jlq — vp) 



dkdqdp 



-*-k I ^t^hm* - m- (27) 

Consider the operator (d/dX)^ 1 , which gives the an- 
tiderivativc of the function it works on. Then we have 



-ikX 



i ( d 



/, k \dX J dX ( 
and Eq.J23 becomes 



-ikX 



\dX 



-ikX 



(28) 



qW(q,p)exp[— ik{X — jlq — vp)] 

L(±Y 1 



dkdqdp 



du \dX J 



w{X,^V) (29) 



Using the same simple operations we obtain the rules of 
Ea. (|26|) terms transformation, which we formally desig- 
nate as '— >': 



A. Evolution equations 

Let us discuss the evolution equation for tomogram 
w. Begin with the most general evolution equation for 
density matrix: 



(25) 



Here and throughout the paper we omit the dependence 
on time t, but imply that all functions, describing the 
state (density matrix, Wigner function, tomogram) de- 
pend on time as parameter. We consider the Hamilto- 
nian H = J^iPi /(^ m i) + V(q). To derive the evolution 
equation for tomograms we consider the Moyal evolution 
equation for Wigner function [Tfill45ll4^ |: 



dW 
~~dt 



"<'+5B>-"<'-5S> 



p dW 
m dq 

w = o, 



(26) 



qW(q,p) - ~ (J^ ) t r(X./7.P) (80) 



dW(q,p) ^8 
<9(f > ^~dX~ W 



(31) 



pW(q,p) - 

dW(q,p) 
dp 



d ( d 



dv \8X 
d 



w{X,jl,v) (32) 



(33) 



Successive application of rules (|30I33|) allows us to trans- 
form all powers of q, p and corresponding derivatives in 
Eq. (|26(l . As a result we obtain the evolution equation for 
one random variable quantum tomogram w: 



dw ft dw . 
dt m dv 





1 _ d 
— v 

2 dX 



(34) 
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B. Stationary states 

For the stationary state with definite energy we can 
turn from the time-dependent Schrodingcr equation 125|) 
to the eigenvalue equation: 



Hp E = PeH = Ep E 



(35) 



Applying the transform (|llfl we obtain the following 
rules of transition from the equation for density matrix 
to equation for Wigner function: 



d 2 p{q,q l ) 
dq 2 

V{q)p(q^ 



After that, using (|30I33() . we have the eigenvalue equation 
for the tomogram w with one random variable: 



Nd 

E 



i d 2 



Inij dv 2 \ dX 



d 



1 



a 1 



-fj, 



&77lj ° dX< 







I i _ d d 

ReV \2 v dx~dn\dx 



Ew 



Nd 

E 



dw 



2m i di 



i _ d d ( d 
2"dX ~ dfl \dX 



w 



(37) 



case we consider 'imaginary time' (3 — 1/T (measuring 
T in units of energy). (3 enters as a parameter in the 
density matrix, which is now defined by the equation 



dp(<f,q",P) 
d(3 



H s ,p(q',q",P), 



(40) 



where index <f in H^i shows that the Hamiltonian acts 
only on those variables. 

Now the transition to the tomograms w\ or w is 
straightforward. We just use the same rules, as in 
the derivation of evolution equation l|34|l and eigenvalue 
equation l|37|l . Then the evolution equation in imaginary 
time (3 for w is given by 



(36) _dw = 




1 



U 2 



-p, 



8m,- ^ dX 2 



w + 



ReV —— - — 



Nd 



E 



Pj dw 



i 9 
ImV 2 V dX 



d_ 

dp 



d 

dX 



(41) 



Initial condition for Ea. ll4l7|) is p(cf,q",/3 = 0) = 
S(q > — q"). It corresponds to constant Wigner function 
(see Ea. ljTT|) '). Using Ea. (fT5|) we see that the tomogram 
w for (3 = must have the delta-function form, equal 
zero everywhere, besides the point p, v = and constant 
in X direction in that point. 



C. Quantum transitions 

In general, there is a possibility of transition between 
the quantum states. Consider two states, designate them 
a and b. The probability of transition from state a to 
state b is P ab = Tr(p a p b ) = J p a (q l ,q")p b (q" ^q^dq 1 ' dq" '. 
In terms of the Wigner formalism this can be rewritten 
as 



P ab = (2n) Nd J F w ^{q,p)F w ^{q lP )dqdp, 



(38) 



and recalling the connection of Wigner function with to- 
mograms w\ and w Ijltjl I14|) we easily get the following 
expressions for P ab in tomography approach: 



w (X, p, v)w (Y, -p, -v)e t(A + 1 >- 



(27T) 



Nd 



(39) 



D. Tomographic map in temperature-dependent 
processes 

The tomographic representation can be analogously in- 
troduced for the systems with temperature T ^ 0. In this 



E. Average values calculation 

Developing the 'center-of-mass' tomography formalism 
we must provide the rules of average values calculation 
to complete the picture. Using the density matrix to 
describe the state of the system we can obtain the average 
value of some operator A as 



(A)^Tr(pA), 



(42) 



where we choose Tr(p) = 1. 

Here it is again convenient to begin with the Wigner- 
Moyal formulation of quantum mechanics. In its frame- 
work to calculate the average value one deals with the 
Weyl symbol A w (q,p) H3 of operator Afap) (see |48ll49l | 
for review): 



(A) = / A w (q,p)W(q,p)dqdp, 



where the Weyl symbol is given by 



A W (q,p) = I Tr{A{q,p)e* q+lrip )e 



iZq+iiJP\„-iZq-iVP 



(43) 



d£,drf 
(2n) 2Nd 

(44) 



6 



Expression for the average values in one random vari- 
able tomography formulation is obtained using the con- 
nection between w and Wigner function l)14[l: 



(A) = 
A(ft, v) 



e 1 w{X,fi,v)A(p,v)dXdp,dv, 



A w (q,p)e 



dqdp 
(2w) 2Nd 



(45) 
(46) 



If considered operator depends on coordinates q or mo- 
menta p only, Weyl symbols have the same form as cor- 
responding operators in coordinate or momentum repre- 
sentation. Operator A(q) is A(x) in cc-coordinate rep- 
resentation, then its Weyl symbol A w (q,p) is equal to 
A(q). The same is valid for momenta-dependent oper- 
ator: B{p) is B(y) in y-momentum representation, and 
B W (q,p) = B(p). 

Consider an operator A(q), depending on coordinates 
only. For momenta-dependent operators all equations are 
the same, provided /i is replaced by v, and vice versa, be- 
cause the pairs q, [2 and p, v enter the equations connect- 
ing the tomogram w with Wigner function symmetrically. 
Integration over z7 in Eg. 1)45(1 for operator A(q) gives the 
delta- function 5(y). Then we have: 

(A) = J A w {q)e-^-x) w{x ^^ =Q) d _^^ m 

It is often necessary to operate with the one-particle and 
one-dimension operators. Then, quite generally, we can 
consider an operator A(qi). Corresponding average value 
is given by 



(A) 



A w (X)w(X,fj ll = l,fl=0,0)dX, (48) 



where fl designates all [ij except the specified [L\ . 



IV. EXAMPLES 

In this section we introduce several examples of to- 
mographic map for many-particles quantum states. For 
simplicity, here we do not regard symmetry over parti- 
cles exchange. Permutations properties are considered in 
Sec. El 



A. Gaussian states 

Quite simple is the case of pure state and wave function 
of Gaussian form. This can be the ground state of the 
system of independent oscillators, as well as coherent or 
squeezed states, or any many-dimensional Gaussian wave 
packet. Such wave packet can be created due to paramet- 
ric excitation of multimode vacuum state of electromag- 
netic field |5lL e.g., in the framework of nonstationary 
Casimir effect 531. 



Consider the pure state with the wave function ^(q) 



lifA^iilj): where 



■4>j{q) = {A J /n) 1 ^e-^-^- x ^ 2+ty ^ 



(49) 



The only mathematical fact we need here is that the 
Fourier transform of a Gaussian is Gaussian. Then, using 
Eq. we immediately obtain the Wigner function as a 
product of Wj(qj,Pj), where 

Wjiq.p) = e-^ii-^Y' '■ e -B j {p-y j y \A j B j f/ 2 /tt, (50) 

and for states 1|49[1 Bj = 1/Aj. For the set of parameters 
x, y, A, B applying Fourier transformation (jT3|l to lj5U|) 
we have: 

- Gauss {X^,u) 



w 



-(X-px-uy) 2 /C 



(51) 



where C = E; = i(^/^ +^/^). 

Thermal density matrix of independent oscillators is 
also Gaussian, but it is not a product of wave functions, 
as the state is not pure. Still it is a product of density 
matrices of individual oscillators (see, e.g. |53(): 



_ j 2A ] {B, J - 1) A 1 fB,( Q 2 +^)-2ga'l 



(52) 



where Aj — mujj /(2sh(uij(3)) and B — ch(ujj(3). Omit- 
ting the straightforward calculations, we obtain the to- 
mogram w in the following form: 



w^{X,^V) 

Nd / 



-X 2 /D 



D 



E 

i=i 



2A,{B 3 



2v]A 3 {B 3 




B. Fock states 

The Fock states of light (the eigenstates in representa- 
tion of photons number) correspond to ground or excited 
states of multimode oscillator. The state is labeled by 
vector n of integer numbers and wave function has the 
form: 



rV4 



m A 



(55) 



where H m is the Hermit polynomial of rn-th order. To 
obtain the tomogram for such state we use the following 
facts. First, coherent state of an oscillator is described 
by the Gaussian wave function and, correspondingly, by 
the Gaussian tomogram (see Ea. (l5T)> ). Coherent state is 
labeled by complex vector a = a + ib and parameters 
of Gaussian wave function in coordinate representation 
(|49|l are xj = \[2aj and y 3 — —\/2bj. Second, the wave 
function of coherent state (for simplicity, one dimension 
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is considered here) is expanded in the basis of Fock states 
as 



\a\ 2 /2j2 



a" 



71=0 



(56) 



which is connected with the expression for generating 
function of Hermit polynomials: 



-oi+2 



n=0 



(57) 



Expanding the tomogram of coherent state in Hermit 
polynomials and wave function of coherent state in cor- 
responding integral expression in wave functions of Fock 
states we have 



w n (X,fl,if) 



Nd 

n 

3=1 



. Nd 



3 



-dX 



(58) 



For example, for N = 2, d — 1 and states with 
nx,n 2 equal to or 1 (denoted (nx,n-i)) the tomograms 
w(X, fix, p2, v\, v<i) have the forms 



(o 
w 2 



„(o,i 



o) exp [-X 2 /C] 

CT \2C 2 X 2 + C X C 2 + C?)e- S 



ttCi 



C 5 /2 



(1>1) _ 4^^e-~ 



2 

X 4 



X 2 C\ + Cf - 4Ci(72 3 
_ + 4 



(59) 
(60) 

(61) 



C 2 C dC 2 
where C\ = /if + v\ , C 2 = ^ + ^ f and C = Ci + C 2 



V. SYMMETRY PROPERTIES WITH RESPECT 
TO PARTICLES PERMUTATIONS 

Consideration of identical particles exchange imposes 
the restrictions concerning the possible form of the state- 
describing functions. In this section we discuss the corre- 
sponding properties of one-random- variable tomogarphic 
map (see [54j for permutation symmetry properties of the 
symplectic tomogram). 

Further we use the following notations. A vector with- 
out index a has Nd components, vector with index dj 
denotes the set of some values, corresponding toj-th par- 
ticle, and consists of d components. A vector a denotes 



the collection of all components of a, except those that 
are specified in the same expression. For example, q in 
the expression ip(qj,q) is the vector of all the coordinates, 
except the coordinates of the j-th particle. 

For particles obeying Fermi or Bose statistics, we have 
the following symmetry properties concerning their per- 
mutations: 



p(qj,qi,q';qi,qj,q") = p(qi,qj,<?;qj,qi,<?') = 

where the upper sign ('+') is for Bose systems, and lower 
sign ('— ') is for Fermi systems. Note that 'entire' par- 
ticles permutation (two particles exchange both q and 
q' variables) corresponds to sign conservation for both 
Fermi and Bose statistics: 

P % c, i') = ptfi, in (63) 

In the expressions for obtaining the Wigner function 
form density matrix and tomogram w from Wigner 
function l|13fl we can exchange the integration variables 
(Hj <-* Uj, etc.), then we immediately have: 

W(q*j,q\,q*;pj,pi,p) =W(qi,qj,q;pi,Pj,p), (64) 
w(X; j2j,pi,fl; Pj,Vi, v) = w(X; pi, pj,p; v^v^V) (65) 

We see that there is no distinction between Fermi and 
Bose statistics when the particles exchange 'entirely', i.e. 
q and q' in the density matrix, q and p in Wigner function 
or in i» are permuted simultaneously. The distinc- 
tion appears when not all the variables, corresponding 
to the considered particles, are permuted. When we use 
the density matrix, Fermi and Bose statistics differ only 
in the sign ±1, which appears after the permutation of 
either q\, or qV, qt'. For the Wigner function and tomo- 
gram this difference is expressed in far more complicated 
manner, through the integral transforms (see correspond- 
ing formulae for the symplectic tomography in |54[). 

First, regard the permutation of ql, qj or Pi,Pj for the 
Wigner function. Again exchanging the integration vari- 
ables in Ijlll) we come to 

W{qj,qi,q;pi,pj,f) = W Pj,Pi,f) (66) 

The same considerations lead us to the similar expres- 
sion for W. 

w(X; flj,Ui,u; V tl Vj, v) = w(X; p h p,j,p; v s , u i} v) (67) 

Then it is enough to develop the formulae for coordi- 
nate (Wigner function) or p, (tomogram) permutations 
only. Corresponding integral expressions has the follow- 
ing form: 
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W(qj,qi,q;pi,pj,p) = I K w (xi, xj, y it yj, q it qj,Pi,Pj)W(xi, Xj, q; y it yj^dxidxjdyidyj, (68) 
w(X, m, vu Vj,v) = j K{X, jli, Vj ; Y, ffi, rJj)w(Y, rji, ff j ,i/)dYdiidi j dffidff j , (69) 

and kernels are given by 

K w {x u x h y u y h qu QjJuPj) = ± (J^j 6{2i + xj - £ - q 3 )S(ti + y, - p, - pf J .) e *[(«-«f)Wi-ft)+(^-^)(ft-?(¥o) 
K(X, jli, flj, Vi, Uj;Y, ii, ij,f]i, rjj) = 

± J J^TT^ + ij ^ WW + % ^) e -^ fc ^- y )- fe2 / 4 [^-^)W--*)+te-^)(^-^)]>dfc (71) 



VI. CONCLUSION 

We studied in details the version of tomographic map 
of the density matrix and Wigner function for which the 
quantum state of multimodc system is associated with 
probability distribution function. This function depends 
on one random variable X and 2iV<i real parameters (real 
iVd-vectors and v) and it determines the quantum 
state completely. It means that provided this probability 
distribution function is known one can reconstruct the 
Wigner function of the system state and corresponding 
density operator. The random variable X can be inter- 
preted as the system " center of mass" coordinate consid- 
ered in specifically rotated and scaled reference frame in 
the complete phase space of the system. Real parame- 
ters (vectors jl and v) determine this rotated and scaled 
reference frame. 

Information contained in the introduced tomogram (w) 
is the same as that contained in the symplectic tomo- 



I 

gram (wi), which depends on larger number of variables. 
It corresponds to the fact that the tomograms have high 
symmetry properties. By means of the symmetry opera- 
tions one can reconstruct the dependence of the function 
on larger number of variables starting from initial func- 
tion with smaller number of variables. 

We have constructed the quantum evolution equations 
and energy level equations for the introduced "center of 
mass" tomogram. Example of multimodc oscillator and 
symmetry properties of the tomogram for identical par- 
ticles (fermions and bosons) were discussed in details. 
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